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Abstract 

A discrete-dynamics model, which is specified solely in terms of the sys- 
tem's equilibrium structure, is defined for the density correlators of a simple 
fluid. This model yields results for the evolution of glassy dynamics which 
are identical with the ones obtained from the mode-coupling theory for ideal 
liquid-glass transitions. The decay of density fluctuations outside the tran- 
sient regime is shown to be given by a superposition of Debye processes. The 
concept of structural relaxation is given a precise meaning. It is proven that 
the long-time part of the mode-coupling-theory solutions is structural relax- 
ation, while the transient motion merely determines an overall time scale for 
the glassy dynamics. 

PACS numbers: 64.70.Pf, 61.20.Lc 
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I. INTRODUCTION 



Glass-forming liquids exhibit a dynamics which appears anomalous in comparison to the 
one of conventional condensed matter. The characteristic time scale r for the motion can be 
several orders of magnitude larger than the natural time scale of, say, normal liquid dynamics. 
The scale r is extremely sensitive to changes of control parameters such as temperature T or 
density n; for example, a change of T by 10 degrees may imply a change of r by a factor 100. 
Furthermore, decay of correlations or spectra are stretched over dynamical windows of time 
t or frequency to, respectively, of many decades in size PJ. The evolution of the anomalous 
dynamics upon cooling or compressing a liquid has first been studied comprehensively by Li 
et al. for the mixed salt Cao.4Ko.6(NOs)i.4 (CKN) and by van Megen and Underwood 
for a hard-sphere colloid 0. The former determined depolarized-light-scattering spectra 
for a four-decade frequency window and the latter measured photon-correlation curves for 
an eight-decade time window. In recent years the mode-coupling theory (MCT) for the 
glass transition has been developed [[|. It deals with a mathematically well-defined model 
for anomalous dynamics, which results from a bifurcation point; its novel features are due 
to the interplay of non-linearities with divergent retardation times. Results of this theory 
have been tested against experiments, for example in Refs. The outcome of these 

tests qualifies MCT as a candidate for a theory of the anomalous dynamics in glass-forming 
liquids. 

The equilibrium structure of a classical system is determined by the ratio of interaction 
potentials and thermal energy via the Boltzmann factors. It is independent of the par- 
ticle masses or other inertia parameters. It is the same for a conventional system, ruled 
by Newton's equations of motion, and for a colloid, whose time evolution is controlled by 
a Brownian dynamics. The simplest information on structure is provided by the struc- 
ture factor S q oc (|f?,f| 2 ), where g^ are density fluctuations of wave vector q; here and in 
the following q — \q\ denotes the vector modulus, and () abbreviates canonical averaging. 
The indicated anomalous dynamics in glassy systems is not related to anomalies of the 
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equilibrium structure; S q is a smooth function of n, T and q throughout the whole liquid 
regime. The simplest information on structure dynamics is provided by the density correla- 
tors $ q (t) = (Qq(t)*0q) / (\Qq\ 2 ) ■ The evaluation of these functions is the main theme of MCT. 
In this paper we will restrict ourselves to simple liquids and to the idealized version of the 
MCT [§. 

Two propositions shall be considered. First, the anomalous dynamics is solely determined 
by the equilibrium structure, i.e., by the potential landscape in the configuration space. 
Second, the anomalous dynamics can be described by a superposition of Debye-relaxation 
processes: 

$ ?(*) = I]PgJ eX P(-7g,^) , Pq,j > , lq,j > . (1) 
3 

These propositions provide a precise meaning to the statement that the anomalous dynamics 
is structure relaxation. 

Density correlators in conventional liquids exhibit the short-time expansion $ g (i) = 
1 — \{£l q i) 2 + 0(t 3 ), where Q 2 = v 2 q 2 /S q , with v 2 = (ksT/m) denotes the bare phonon 
frequency |J. This expansion contradicts Eq. ([1]). Furthermore, the thermal velocity v 
depends on the particle mass m and is therefore not an equilibrium quantity. It is necessary 
to formulate the propositions more precisely to avoid a conflict with the cited short-time 
behavior. Thereby some insight in the physics of glassy dynamics shall be provided. 

II. BASIC EQUATIONS 

Within the Zwanzig-Mori formalism one can derive the equation 

d 2 <S> q (t) + {l 2 q $ q (t) + f M q (t - t')d t ^ q {t')dt' = , (2) 
J 

where M q (t) is a correlation function of fluctuating forces 0. The kernel M q (t) can be 
split into a regular part M^ cg (t), dealing with conventional liquid-state dynamics, and a 
contribution Q 2 m g (t), dealing with slowly fluctuating forces due to sluggishly moving struc- 
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ture: M q (t) = M r q eg (t) + Q 2 m q (t). Normal-state dynamics would then be obtained from the 
equation V<& q (t) = 0, where the abbreviation is used: 

v%{t) = d 2 t <$> q (t) + n;$ ff (t) + f M^(t - t')d t ,$ g (t')dt' . (3) 

The subtleties of MCT are due to the approximation of m q (t) as mode-coupling functional 
T ■ 

m q (t) = F q ($ k (t)) = ^ kp V q ^ k {t)%{t) . (4) 

Here the positive coefficients V q ^ P are given in terms of S q , i.e. by the equilibrium structure. 
The MCT equations of motion are 

T>%(t) + n 2 g f m q {t-t')d t ^ q {t')dt' = . (5) 

The theory is specified by S q and M r q cg {t) where neither quantity reflects glass-dynamics 
anomalies. The Eqs. (|3|)~(iD are regular and fl q , V q ^ P , M q eg (t) depend smoothly on control 
parameters. For mathematical convenience the wave-vector moduli are discretisized to a 
set of M values up to some cutoff g max - Thereby MCT deals with M nonlinear integro- 
differential equations of the Volterra type, which are coupled via the functional T q . A 
review of the derivation of the MCT equations, in particular of the explicit form of V q ^ p , 
and citations of the original papers can be found in Ref. 0. 

Let us consider the hard-sphere system (HSS) as main example for the following demon- 
stration of our results. Its equilibrium state is controlled by the packing fraction (p = 7md 3 / 6. 
As unit of length the particle diameter d is chosen. Wave vectors are considered up to the 
cutoff g max = 40 and we use M = 100 equally spaced q values. The structure factor is 
evaluated within the Percus-Yevick theory 0. Two models for the regular dynamics shall 
be studied. In the first one the regular kernel is dropped, so that 

V^%(t)=d^ q (t)+n^ q (t) . (6a) 

This is a model with Newtonian dynamics where the transient deals with oscillations. The 
second model is obtained by dropping inertia effects and replacing the regular kernel by a 
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g-independent friction term v. This is a model for a colloid where the transient deals with 
relaxators 

p (2) $ 9 (t) = vd t %{t) + n 2 q $ q (t) . (6b) 

Another set of examples shall be formulated for a so-called schematic model (SM), 
dealing with a single correlator $(t) only [§]. The mode-coupling functional is specified as 
= t>i$(t) + t> 2 $(t) 2 , where t> i > , v 2 > denote coupling constants. Despite its 
apparent simplicity, this M — 1 model reproduces some generic features of the MCT so 
faithfully, that it has been used as basis for a quantitative description of the evolution of 
anomalous dynamics of glycerol within the full gigahertz band 0. For the regular motion 
we again consider the oscillator model, Eq. fl6"aj), and the relaxator model, Eq. (|6b"|) . In 



addition we also study a soft-mode contribution to the regular dynamics. Here the transient, 
described by T> = T>^\ is given by a regular kernel in Eq. (||), which obeys the oscillator 
equation 

(d? + v d t + n 2 )M^(t) = 0, (7) 

with initial conditions M reg (0) = Co , d t M reg (0) = 0. Choosing a low soft-mode frequency Qq 
we shall investigate the interference of slow regular dynamics with anomalous dynamics. 

III. ANOMALOUS MCT DYNAMICS 

The solutions & q {t) of the MCT equations of motion decay to zero for large times as 
expected for a liquid, provided control parameters like n or 1/T are smaller than some critical 
value n c or 1/T C . However, for n > n c or 1/T > 1/T C the solutions describe ideal glass states 
characterized by a positive Debye- Waller factor f q = — > oo) 0. For the HSS one finds 
the bifurcation point for (p c w 0.52 [K|, which is not too far from the experimental value 
0.58 [Q. The evolution of the HSS dynamics for wave vector q = 10.6, calculated for the 
oscillation transient (|5"aD, is shown in Fig. 1; and Figs. 2 and 3 exhibit the corresponding 
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fluctuation spectra &g(u)) and susceptibility spectra Xqi^) — ^Qqiuj), respectively. Here 
<&q(u)) = /o°° cos((jt)$> q (t)dt is the Fourier cosine transform of the correlator. The wave 
vector q = 10.6 is located close to the first minimum of the structure factor. For the SM 
the weak-coupling liquid regime is separated from the strong-coupling glass regime in the 
v\ — t>2-plane by the parabola v\ = (2A — 1)/A 2 , v\ = 1/A 2 , 1/2 < A < 1 |. Figures 
4-6 exhibit the evolution of the dynamics for the soft-mode model, Eq. (0), upon crossing 
the transition line at A = 0.7. The Figs. 1-6 exhibit slow, control-parameter-sensitive, 
stretched dynamics. Obviously, the results for the HSS are similar to the corresponding 
ones for the SM and this exemplifies the relevance of the latter for a discussion of MCT 
findings. The shown anomalous dynamics can be understood from the asymptotic solutions 
of the MCT equations for long times and low frequencies near the transition, as explained 



in the preceding literature, e.g. in Ref. \'.\ I 



There is a qualitative difference for the short-time dynamics of the HSS between the 
results shown in Fig. 1 and the ones shown in Ref. [|TTJ for the dynamics calculated for the 
relaxation transient, Eq. (|6l5|). The correlators of the second model decrease monotonously 
while the correlators in Fig. 1 exhibit oscillations. These get more and more damped if (p 
increases towards tp c . Increasing ip above ip c , the oscillations become again rather pronounced 
in the glass state. These yield the oscillation bumps for the spectra &q{uj) for u > 10 in Fig. 
2, while the spectra for the colloid model decrease monotonously with increasing frequency 



11]. Similar differences occur for the SM. The oscillations in Fig. 4 are different from the 



ones calculated in Ref. JT2[ for the transient (6a); and the correlators for the relaxation 
transient, Eq. flSb]), do not show oscillations at all 0. The fluctuation spectra in Fig. 5 
exhibit two oscillation bumps for u > 0.1, the spectra calculated with Eq. (|6a]) have one 
bump for the glass states ||12|| , and the ones for the relaxation transient ( |6b| ) have no bumps 
at all 0. However, the long-time parts of the decay curves referring to the same mode- 
coupling functional T q coincide for the various models for the transient. This holds provided 
an appropriate overall shift of the curves parallel to the logarithmic abscissa is done in the 
figures. The statement is demonstrated in Figs. 7, 8 for the two HSS models and in Figs. 



9, 10 for the three mentioned SMs. 

The findings in Figs. 7-10 suggest the following precise formulation of the first proposi- 
tion. The transient dynamics determines a time scale to so that for t ^> to the correlators 
can be written as 

$ 9 {t) = F q (t/t ) . (8) 

Here F q is independent of the transient as quantified in Eq. (|]) by Q q and M q eg (t). The long- 
time decay, its sensitive dependence on control parameters and also its dependence on q, is 
given by the master functions F q , which are determined by the mode-coupling functional 
T q . Of course, the transient may also depend on control parameters, and this leads to a 
smooth variation of to- This effect results in the parallel shift of the n = 2, e < 0-curves in 
Fig. 7 and similar offsets of the n = 1-curves in Figs. 9, 10. The proposition is meant as 
asymptotic result for the approach towards the critical point so that to is the limit result at 
the bifurcation singularity. 



IV. THE REGULAR PART t OF THE RELAXATION SCALES 

The first step of the proof of the first proposition is based on the leading and next- 
to-leading long-time expansion of the critical correlators, i.e., of the MCT solutions at the 
bifurcation point: 

= n + h q (t /t) a {l + [K q + K(a)}(t /t) a } . (9) 

Here terms of order (to/t) 3a have been dropped. There are straightforward formulas for the 
evaluation of the critical form factor f q > 0, the critical amplitude h q > 0, the correction 
amplitude [K q + n(a)], and the critical exponent a, < a < 0.5, from the mode-coupling 
functional ||11|| . These parameters are equilibrium quantities. For the HSS one gets a = 



0.312, /fo g = 0.417, h 10S = 0.642, K ws + k(o) = -0.185. For the SM one finds a = 
0.327, f c = 0.3, h = 0.7, K = 0, k(o) = 0.0528. The transient, no matter how complicated 
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M q eg (t) in Eq. (|j) may be, merely enters via the scale i - Specializing Eq. (§) to the 
critical point, one identifies to with the scale in the proposition. There are two diverging 
time scales hidden in <& q (t), which govern its sensitive control-parameter dependence ||. 
Formula ([8]) implies, that to enters these scales as a prefactor. The singular part of the 
scales is determined by the master functions F q (t). 

The found results are demonstrated in Figs. 11 and 12. For all models of Sec. |T|, Eq. 
(g) was matched to the numerical solution at the critical point. This determines the values 
to, cited in the figure captions. Notice in Fig. 12 that a larger percentage of the decay of 
is described by the law @ for the relaxator model than for the oscillator models. Slow 
transient oscillations can destroy the short-time part of the critical fractal decay. Given an 
upper time cutoff t max , the transient oscillations can be modeled so that they destroy the 
critical decay for t < t max completely. However, the decay for t > t max is robust. One can 
construct models so that there is a fully developed a process, i.e., a stretched decay of 3> g (i) 
from f° to zero, without a critical decay precursor. 

V. A DISCRETE DYNAMICS MODEL 

The second step of the proof of Eq. (|j) is based on the Fourier-Laplace transform 
of the equations of motion (|), <g): $ q (u) = -l/[cu - 9?J[uj + M r q ^{uj) + tt 2 q m q (cu)]]. 
Here we use the convention for the transform of some function G(t) to G(u): G(uj) = 
i Jq°° exp(izt)G(t)dt, z = u + iO. At the transition point one derives from Eq. @ a divergent 
small-frequency mode-coupling kernel m q {uj) = [—Fg^f^/u] + 0(l/u 1 ^ a ). Because of con- 
tinuity one concludes that in the asymptotic limit of small frequencies and small separations 
from the critical point, the regular function u + M v q eg {uo) can be neglected in comparison to 
VL 2 q m q {uj). Hence, in this limit the correlators obey m q (u) — & q (uj) = Lum q (uj)<& q (uj) 0. Let 
us assume that can be continued as a function F q (t) to small times so, that the formu- 
lated equation holds for all times and frequencies: [N q (ui) — F q {u)]/uj = N q (u)F q (uj) , N q (t) = 
J~q{Fh{t))- The details of the continuation are of no concern, since we are not interested 
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in the short-time transient. Backtransformation yields the set of M implicit functional 
equations for the M functions F q (t): 



J [N q (t') - F q (t')]dt' = jf N q (t - t')F q (t')dt' , (10a) 

N q (t)=F q (F k (t)), q = l,...,M. (10b) 

Equations of a similar form have been studied before in some different context fl3| , |l4l , and 
we shall adopt some of the tricks of the preceding work to deal with the present problem. 
Notice that Eqs. ( |10"D cannot define a time scale, since they are scale invariant. With F q (t) 
also F q (t) = F q (x • t) is a solution for all x > 0, q = 1, . . . , M . 

The integrals in Eq. ( |10a| ) shall be written as Riemann sums formed on a time grid of 
equal step size 5. The sums shall be determined by the values of the functions in the middle 
of the intervals 

gU=F q ((i + 1/2)5), z = 0,l,.... (11) 
The sums can be regrouped so that Eqs. (|T0|) read 

9l m) =l*(jg? ) ,9?\---,9i n) ), (12a) 
where the functional I q is given by 

{m— 1 m— 1 ^ 

(i - 4 0) )^(4 m) ) + E - of) - E ^{gt'^gf /[i + r q lsf>)\ ■ (12b) 
i=0 i=l ) 

Formula (12a) can be considered as an implicit equation to determine in terms of the 
for % preceding m: 

gM=T q (g®g®...Jr 1) )- (13) 

The explicit solution for g^ m \ i.e., the construction of the functional T q , is done by the 

,M 

)q,t i 



following procedure. One defines a sequence of approximands g^f ,£ = 0,1,..., so that 



lim^oo flfjy = g { q m) . Here g^ +x = X q {gf\ . . . , g { ™ X \g { ™z), and the start is chosen as 
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9ko = S'i" 1 " 1 ''- Notice that the scale invariance is reflected by the fact, that the step size 5 
does not occur in Eqs. flT2|) , ([13|) . 

The result ([K|) can be interpreted as an iterated mapping with memory or, because of 
Eq. fllTD , as the definition of a discrete dynamics with retardation. The M ■ m numbers 
g\, ,k = 1, . . . , M, % = 0, . . . , m — 1 determine the M numbers g q m \ A sequence of values 

is created, once that the M initial values g^ are specified. The arbitrariness of the scale 
is hidden in the one of the choice of the initial condition. 

By construction one expects that the sequence defines a solution of Eqs. (|10|) . To show 
this explicitly one can read all equations backwards. As an approximand for the F q (t), step 
functions are defined by F & q (t) = gf for i < t/5 < (i + 1). Obviously, lim$_»o = F q (t) 
solves Eqs. (0), and these functions agree with the master functions in Eq. (|8|) for t large 
compared to the time scale to- The limit 5 — > for t > to is trivially taken by using Eq. 
(PI) only in the limit i — >• oo. In practice, a large %q is chosen and the g^' for i < io are 
considered as transient of the discrete mapping. For the prescribed small step size 5, 5i has 
to be smaller than the lower bound of the time window to be studied. The found results are 
then independent of the indicated details, up to an overall time scale t^. To cope with the 
stretching of the dynamics over many decades, we applied a decimation procedure ||14|| . 

The discrete dynamics and thus the master functions F q have been constructed in Eqs. 
(PI), (PI), (|13D solely in terms of the mode-coupling functional T q . Hence the function F q 
in Eq. (f|) is a quantity defined by equilibrium distributions. This finishes the proof of the 
first proposition. Figures 13 and 14 exhibit the construction of the master functions F q (t) 
for the HSS and the SM. The corresponding results are included as dotted lines in Figs. 
7-12. 

Some side remark shall be added. Equations (10), specialized to the functional at the 
critical point and complemented by the initial condition F q (t = 0) = f q , are equivalent to 
the equation for the MCT a-relaxation master function ||. The specified iterated mapping 
used with initial condition gj® = f q — 5*h q , 5* > 0, 5* — * 0, is a very simple and most 
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efficient algorithm for the evaluation of F q (t). 



VI. COMPLETE MONOTONICITY 

Let us combine Eqs. fl5|), (|6l5| ) to the MCT equations of motion for M relaxators, specified 
by the decay times r q = v /Q? q 

T q d t $ q (t) + + f m q (t - t')d t ^ q (t')dt' = . (14) 

Suppose that the kernels m q (t) are given as a superposition of L Debye processes: m q *\t) = 
Z)j=i exp(— r^t), where /£ > and the rates are labeled so that < < Y 2 q < ■ ■ ■ < F q . 
The Fourier-Laplace transform m q *\uS) is a meromorphic function with L simple poles at 
— ir^. Substitution of this result in the Fourier-Laplace transform of Eq. (|i"4|) yields the 
solution as meromorphic function. Elementary discussion brings out that $*„{uS) has 

exactly (L + 1) poles which are located on the imaginary axis at, say, (—ij q ); the poles of 
m*(u) separate those of i.e., > j q ' > T q for j — 1, . . . , L — 1, 7^ > r£, > 

7° > 0. The residues p> of the poles are positive. Hence the solution is a superposition of 
(L + 1) Debye processes = Ej= 4 ex P( _ 7^)- 

Neglecting the kernel m g (t), the solutions of Eq. ( |I4"D are Debye processes = 
exp(— t/r q ). The relaxators are coupled by the mode-coupling functional (4), which is a 
polynomial with positive coefficients V q ^ v . Hence = mW(t) is a sum of a finite 

number, say L, Debye functions. Substitution of this kernel into Eq. ([14]), leads to a solution 
$^(t), which according to the preceding paragraph is a sum of (L + l) Debye contributions. 
This can be used to define a new sum of Debye processes m q 2 \t) = Continuing 
one constructs a sequence of approximands m q n \i), & q n \t), n = 0,1,2,... of the type 
formulated in Eq. ([[]). In Ref. ]TB| it was shown: the sequence converges uniformly 



towards the unique solutions of Eq. (|1~4" ) and this solution is completely monotone, i.e., it 
obeys (-d/dt) e $ q (t) > , £ = 0, 1, . . . . 

According to Bernstein's theorem [I5|, every completely monotone function can be rep- 
resented as Stieltjes integral 
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%(t) = / exp(-7t)da 3 ( 7 ) , (15) 
J o 

where the measure 015(7) is an increasing function of the rate 7. The theory of these integrals 
implies the following. For every finite interval of positive times t, 0<ti<t<t 2 <oo and 
every positive error margin 77, one can find a set of numbers (A > , 7^ > ,j = 1, . . . , N, 
so that 

l^(*)-Ef =1 4 ex P(-T^)l<^- (16) 

The solutions of Eq. flTj|) can be approximated arbitrarily well by a finite sum of Debye 
processes. 

Because of Eq. (^j) the solutions of any MCT model are given for t ^> to by the transient 
independent function F q (t/to). This holds in particular for the solution of Eq. (|Hp. Choos- 
ing ti larger than t one can replace $,(*) by F q (t/t ) in Eq. flT|). Thereby one obtains the 
precise formulation of the second proposition and its proof. 

VII. CONCLUSIONS 

The dynamics of a classical system deals with the orbits in phase space. Collision events 
and vibrations are the elementary bits building the motion on microscopic scales. In dense 
liquids there is the cage effect |I7| : the system gets trapped in phase space pockets for 
long times. This implies a separation of a low-frequency contribution to the spectra from 
the normal-motion band. This contribution is the anomalous dynamics dealing with the 
motion from one pocket to the other. The first proposition leads to the following picture. 
The anomalous dynamics reflects the statistics of orbits in configuration space. After coarse 
graining of the time over intervals of microscopic size, the normal condensed-matter dy- 
namics does not play a role anymore; it merely sets the scale to f° r the exploration of the 
potential landscape in the high-dimensional configuration space. It does not matter, for 
example, whether this exploration is done as prescribed by Newton's equations of motion 
or by Brownian dynamics. Consequently, the correlation functions for the coarse-grained 
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dynamics are governed, up to an overall time scale to, only by Boltzmann factors, i.e. by 
equilibrium distribution functions. 



The statistics of orbits is studied in the theory of generalized Brownian motion [JT8 
From this theory one expects as generic results time fractals like the one formulated in 
Eq. ([5|). These fractals result from the mapping of the orbits on the time axis as achieved 
by correlation functions. They reflect nontrivial Hausdorff dimensionalities of sets defined 
for return-time and waiting-time distributions. MCT should be viewed as a mathematical 
model allowing an explicit evaluation of such orbit statistics via the master functions F q in 
Eq. (H) or via the discrete-dynamics procedure of Sec. |V]. 

Debye's relaxation law is the paradigm for a dynamics without memory. Conventionally, 
it is derived by assuming random forces with a white-noise spectrum [§,[17]] • The derived 
Eqs. flU), (|16D formulate a similar picture for the anomalous dynamics in glass-forming 
systems. Memory effects are irrelevant for the coarse-grained orbits through the potential 
landscape. At the first glance this finding appears as a contradiction to MCT, which deals 
with a retarded fluctuating force via the integral term in the equation of motion Indeed, 
it is this memory term, which renders the MCT bifurcation scenario so different from what 
one knows for bifurcations of conventional dynamical systems. The solution of the paradox 
was explained in the first paragraph of Sec. |VI|. A white-noise spectrum is sufficient to 
produce the Debye law, but it is not necessary. If the forces are superpositions of L Debye 
laws, their retarded influence leads to a response, which is also a superposition of Debye laws 
albeit of (L+l) terms. The MCT memory effects do not destroy the complete monotonicity, 
i.e., the possibility for a representation as superposition of Debye processes. But the memory 
effects lead necessarily to a distribution of the relaxation rates, i.e. to relaxation stretching. 

Our result fl3|) implies that the susceptibility spectra can be written as x'qd^) — Xq(ut ). 
Here the master spectrum x 9 (u;) is given by the equilibrium structure and it can be evaluated 
from the Fourier cosine transform F"(u) of the master functions F q (t), discussed in Sec. [V]: 
X q (uj) = (jjF"(oj). This holds, provided the frequencies are sufficiently small compared 
to the characteristic scale 1/to, which is determined by the transient. This conclusion is 
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demonstrated in Fig. |15| for the two HSS models. For n > 6, i.e. for |e| = \(<p — <p c )/<p c \ < 
10~ 2 , and oot < 10~ 2 ' 5 the spectra for the relaxator model (full lines) agree with those for 
the oscillator model (dashed lines), and both agree with the results obtained for the discrete- 
dynamics model (dotted lines). Two features of our findings should be emphasized. Firstly, 
the full lines represent spectra for the relaxator model, but only the part for log 10 ut < —0.5 
is structural relaxation as given by x q (ujto). Due to mode-coupling effects there are non- 
Debye relaxation spectra, which deal with crossover phenomena from transient to structural- 
relaxation dynamics. The corresponding crossover window is much larger for the oscillator 
model; it extends over the window —2.5 < log 10 co>to < —0.5. Here the dashed curves in Fig. 
15| can be described reasonably by an effective power law x'gi ^) = h^(ujto) acfi , where the 



effective amplitude and the effective exponent a e g depend on the control parameter if. 
Secondly, the range of the applicability of Eq. (|8]) can be extended by incorporating the 



smooth control-parameter dependence of to, mentioned at the end of Sec. pT[ Thereby the 
$(t) versus log 10 (t/t ) diagrams, shown e.g. for e < 0, n — 1 in Fig. |9], get a shift parallel 
to the abscissa so that they coincide with the dotted line there. Similarly, negative shifts 
bring the low-frequency spectra in Fig. for e < 0, n = 4 for the two HSS models on the 
master spectra, and positive shifts do the same for the e > 0, n = 4 results. The important 
point is, that the same shifts do the rescaling for all wave vectors q for the HSS. 

The critical decay law (P) is the germ of all analytical results derived within MCT. It was 



first measured for CKN and that by inelastic neutron scattering ||19|| , and by polarized as well 
as depolarized-light-scattering spectroscopy j2(J . Later this part of the anomalous spectrum 
was also explored for a number of other systems like, for example, orthoterphenyl pl|-p3 



The function / + h/t a with / > , h > and < a < 1 is completely monotone and thus it 
can be written as Eq. ([[]). The l/t a -decay leads to a self-similar spectrum $"(u;) oc c<j 1_a , 
which does not define a time scale. The time scales 1/7^ in Eq. (|1|) merely reflect the 
dynamical window, within which the critical law was studied. The representation of the 
critical spectrum by Debye peaks []22| is not an alternative, let alone a phenomenological 



theory, for the original formulation of the discovery HT9^pT| , rather it is a reformulation in 

14 



a manner suggested by Bernstein's theorem. We have shown in Sec. |V| that not only the 
critical spectrum, but the entire low-frequency spectrum can be represented as sum of Debye 
contributions. 
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FIGURES 



FIG. 1. 

Density correlators for the HSS for wave vector q = 10.6 calculated for the oscillator tran- 
sient, Eq. (pa]). The time unit is chosen so that v/d = 2.5. Curve c refers to the crit- 
ical packing fraction ip = ip c . The uppermost curve is calculated for ip = 0.60, and the 
other curves refer to (ip — ip c )/Lp c = e = ±l/10 n//3 , n = 0, 1, . . .. The free-oscillator curve 
$ g (t) = cosQ q t, which refers to e = —1, is drawn only for Q q t < 1.98. The arrow marks the 
time t = 0.00944. 

FIG. 2. 

Density fluctuation spectra for the results of Fig. 1. 

FIG. 3. 

Susceptibility spectra for the results of Fig. 1. 

FIG. 4. 

Correlators for the SM with soft-mode transient according to Eq. ([?]) with parameters 
i> = 0.2Q, Q = 0.5f2, Co = 0.2. The time unit is chosen so, that Q = 1. Curve c refers 
to the critical point v{ = (2A - 1)/A 2 ; v% = 1/A 2 , A = 0.7. The others are calculated for 
v i,2 — ^12(1 + e )> t = ±l/4 n , n = 0, 1, . . .. The arrow marks the time to — 0.0649. 

FIG. 5. 

Fluctuation spectra for the results of Fig. 4. 
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FIG. 6. 

Susceptibility spectra for the results of Fig. 4. 

FIG. 7. 

Curves (1) are e < results from Fig. 1 replotted with scale to — 0.00944. Curves (2) are 
the corresponding results for the relaxation transient, Eq. (|6~bD, replotted from Ref. |yj 
with scale to = 0.425. The dotted lines (3) are the discrete-dynamics-model results from 
Fig. 13, replotted with to — 2675. The labels n are set as in Fig. 1. Curves for successive 
values of n are shifted horizontally by two decades in order to avoid overcrowding. 

FIG. 8. 

The analogous results as shown in Fig. 7 but for e > 0. 

FIG. 9. 

Curves (3) are e < results from Fig. 4 replotted with scale t Q Q = 0.0649. Curves (1) 
(and (2)) are the corresponding results calculated for the oscillation transient, Eq. (|6a]) 
(relaxation transient, Eq. (|6~b|)) replotted from Ref. fl2| (Ref. 0) with scales to& = 0.0440 
(t v = 0.150). The dotted lines (4) are discrete-dynamics-model results from Fig. 14, 
replotted with scale t — 26.25. The labels n are set as in Fig. 4. Curves for successive 
values of n are shifted horizontally by two decades in order to avoid overcrowding. 

FIG. 10. 

The analogous results as shown in Fig. 9 but for e > 0. 
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FIG. 11. 

Density correlators of the HSS for q = 10.6 and the critical value (p c for the packing fraction. 
Curve (1) is replotted from Fig. 1 with to = 0.00944; curve (2) refers to the relaxation 



transient, Eq. (|6"b"|), and it is replotted from Ref. |TT[ with to = 0.425. Curve (3) refers 
to the discrete dynamics-model and is replotted from Fig. 13 with t = 2675. Curve (5) 
exhibits the asymptotic expansion law, Eq. (|9|), and the curve (4) is the leading contribution 
toEq. (1): %(t) =f c q +h q {t Q /t) a . 



FIG. 12. 

Decay curves for the SM for the critical point A = 0.7. Curve (1) refers to the oscillation 



model, Eq. (|6a|); it is taken from Ref. [12] and replotted with t^Vt = 0.0440. Curve (2) refers 
to the relaxation model, Eq. (|6"15D ; it is taken from Ref. [JJ and replotted with t^v = 0.1498. 
Curve (3) refers to the soft-mode model; it is replotted from Fig. 4 with t Q = 0.0649. Curve 
(4) refers to the discrete-dynamics model and it is replotted from Fig. 14 with t = 26.25. 
The dashed line (5) exhibits the asymptotic expansion law, Eq. (P). 



FIG. 13. 

Density correlators of the HSS for wave vector q = 10.6 calculated for the discrete-dynamics 
model with initial condition gjP' = 10. The time unit is chosen so that S = 1 in Eq. (fLip . 
The labeling of the curves is done as in Fig. 1. The arrow marks the time to — 267. 



FIG. 14. 

Correlators of the SM for the discrete dynamics with g(°> = 10. The time unit is chosen so 
that 5 = 1 in Eq. (|TTJ) . The labeling of the curves is done as in Fig. 4. The arrow marks 
the time to = 26.2. 
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FIG. 15. 

Susceptibility spectra for the HSS. The labeling of the curves is done as in Fig. [I]. The full 
lines are calculated for the relaxator transient, Eq. (fTEP, the dashed lines for the oscillator 
model (Eq. (pa])), and the dotted results refer to the discrete dynamics model of Sec. [V]. 
The scales t are listed in the caption of Fig. [7|. 
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